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Abstract: We determine the critical equation of state of the three-dimensional 0(4) 
universality class. We first consider the small-field expansion of the effective potential 
(Helmholtz free energy). Then, we apply a systematic approximation scheme based 
on polynomial parametric representations that are valid in the whole critical regime, 
satisfy the correct analytic properties (Griffiths' analyticity), take into account the 
Goldstone singularities at the coexistence curve, and match the small-field expansion 
of the effective potential. From the approximate representations of the equation of 
state, we obtain estimates of several universal amplitude ratios. 
The three-dimensional 0(4) universality class is expected to describe the finite- 
temperature chiral transition of quantum chromo dynamics with two light flavors. 
Within this picture, the 0(4) critical equation of state relates the reduced temper- 
ature, the quark masses, and the condensates around T c in the limit of vanishing 
quark masses. 



Keywords: 



Dimensions 



QCD, Lattice QCD, Thermal Field Theory, Field Theories in Lower 



Contents 



|T|. Introduction |T| 
|2|. The critical equation of state [5] 



2TT| General properties 



The small-field expansion of the effective potential |7| 



|3|. Parametric representations of the equation of state [8| 
|3.1| Polynomial approximation schemes |8] 
|3l Results ITT 



BL Universal amplitude ratios WA 



1. Introduction 

In the theory of critical phenomena continuous phase transitions can be classified into 
universality classes determined only by a few properties characterizing the system, 
such as the space dimensionality, the range of interaction, the number of components 
of the order parameter, and the symmetry. Renormalization-group theory predicts 
that critical exponents and scaling functions are the same for all systems belonging to 
a given universality class. Here we consider the three-dimensional 0(4) universality 
class, which is characterized by a four- component order parameter, 0(4) symmetry, 
and effective short-range interactions. A representative of this universality class is 
the lattice spin model 

n = -j 3 • sj - ^# • ( L1 ) 

<ij> i 

where Si are four-component unit spins and the summation is extended over all 
nearest-neighbor pairs < ij >. We are interested in the critical equation of state 
that relates the magnetization M = -^(X^^i), the temperature T, and a uniform 
external (magnetic) field H near the critical point T = T c and H = 0. The three- 
dimensional 0(4) model is relevant for the finite-temperature behavior of quantum 
chromodynamics (QCD) with two light-quark flavors. Using symmetry arguments, 
it has been argued that, if the finite-temperature transition is continuous, it should 
belong to the same universality class of the three-dimensional 0(4) vector model 

O § & B- 
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The thermodynamics of quarks and gluons described by QCD is characterized 
by a transition from a low-temperature hadronic phase, in which chiral symmetry is 
broken, to a high-temperature phase with deconfined quarks and gluons (quark-gluon 
plasma), in which chiral symmetry is restored. The main features of this transition 
depend crucially on the QCD parameters, such as the number Nf of flavors and 
the quark masses. Our qualitative understanding of the deconfmement transition is 
essentially based on the expected symmetry-breaking pattern that has been discussed 
at length in the light- [l], |^, |||, f|] and in the heavy- || quark regime. For recent reviews 
see, e.g., Refs. H 0, H, H [K| HU- in the limit of infinitely heavy quarks, i.e. in the 



pure SU(3) gauge theory, the order parameter for the transition is the Polyakov loop. 
The corresponding effective theory turns out to be a three-dimensional -^-symmetric 
spin model, which is expected to undergo a first-order transition. With decreasing 
the quark masses, the first-order transition should persist up to a critical surface (a 
line if only two quarks are considered) in the quark-mass phase diagram, where the 
transition becomes continuous and is expected to be Ising-like JTJ. When the quark 
masses are further decreased, the phase transition disappears and one only observes 
an analytic crossover. In the opposite limit, i.e. for vanishing quark masses, we 
expect again a phase transition that is essentially related to the restoring of the chiral 
symmetry. Its nature can be argued by using symmetry arguments. In the chiral 
limit, since the axial U(1)a symmetry is broken by the anomaly, the relevant global 
symmetry group is expected to be SU(Nf)t X SU(Nf)u. At T = the symmetry 
is spontaneously broken to SU (Nf) L+R with Nj — 1 Goldstone particles (pions and 
kaons). With increasing T, QCD is expected to undergo a chiral symmetry-restoring 
transition, where the order parameter is played by the expectation value of the quark 
bilinear Mj = ((t L qRj)- In the case the chiral symmetry-restoring phase transition is 
of first-order type, it persists also for nonvanishing masses, up to a critical surface 
in the quark-mass phase diagram where the transition becomes continuous and is 
expected to be Ising-like ||. Outside this surface, i.e. for larger quark masses, the 
phase transition disappears and we are back in the analytic-crossover region discussed 
above. In the case the phase transition is continuous, its universal critical behavior 
can be described by an effective three-dimensional theory. 1 Moreover, an analytic 
crossover is expected for nonvanishing quark masses, since the quark masses act as 
a external field coupled to the order parameter. 

In the two-flavor case, Nf = 2, the symmetry-breaking pattern of the transition 
is equivalent to 0(4) — > 0(3). According to universality arguments, if two-flavor 
QCD undergoes a continuous transition, its critical behavior at T c should be that 



1 A finite-temperature e?-dimensional quantum field theory is equivalent to a classical model 
denned in (d+l) dimensions with a finite extent in the "temporal" (d+1) dimension. If the transition 
is continuous, near the critical point the critical modes have a correlation length much larger than 
the temporal extension that can therefore be neglected. Thus, the theory that gives the universal 
features of the transition is effectively d-dimensional. 
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of the three-dimensional 0(4) universality class. Since the quark masses act as an 
external (magnetic) field, a quark-mass term smooths out the singularity and instead 
of the transition one observes an analytic crossover, as in the 0(4) model for non- 
vanishing magnetic field. In the small quark-mass regime the universal features of 
this analytic crossover can be still described by the scaling 0(4) theory. However, 
universality arguments do not exclude that the transition is of first order. In this 
case, it would appear even for small nonvanishing masses. A first-order transition 
should be expected if the breaking of the axial U(1)a symmetry is effectively weak 
at T c . Indeed, in the absence of the U(1)a axial anomaly, the symmetry breaking 
pattern is U{2) L x U{2) R — > U(2) L+R . The corresponding effective three-dimensional 
Landau-Ginzburg- Wilson theory [|l|] does not have stable fixed points, and therefore 
one expects a first-order transition. 2 Since a first-order transition is generally robust 
under perturbations, for a sufficiently small breaking of the axial U(1)a symmetry 
at T c the transition should maintain its first-order nature. This issue has been in- 
vestigated on the lattice, see, e.g., Refs. |TB|, [TJ|, |T5|] . The U(1)a symmetry is not 



restored at T c and the transition is apparently of second order. 3 

In the case of three light flavors, Nf = 3, on the basis of the symmetry arguments 
reported above, it has also been argued that the transition is of first order, up to a 
critical surface in the quark-mass phase diagram, where it becomes continuous and 
it is expected to be in the Ising universality class M. This picture has been also 
supported by a recent Monte Carlo simulation flT7| . See, e.g., Ref. || for a discussion 
of the phase diagram for generic values of the quark masses. A similar behavior has 
been put forward for Nf > 3 [Q. 4 

Since many years the finite-temperature behavior of QCD has been investigated 
by numerical Monte Carlo simulations exploiting the lattice formulation of QCD. The 
available results are substantially consistent with the above-reported picture based 



on symmetry and universality, see, e.g., Refs. |1| [19|, |2(| |2T], g^, g3|. However, 
although the evidence for a continuous transition in two-flavor QCD at T c rs 172 
MeV is rather convincing (see, e.g., Refs. [18, |2~l"ll), conclusive evidence in favor of an 



0(4) scaling behavior in the continuum limit has not been achieved yet. 

An accurate determination of the universal features of three-dimensional 0(4) 
systems, such as the critical equation of state, is important to achieve an unam- 



2 This picture was originally argued on the basis of a first-order perturbative calculation in 
the framework of the e = 4 — d expansion [J. Recently, a six-loop computation within a three- 
dimensional perturbative scheme [l2| put this result on a firmer ground. 

3 The U(1)a symmetry appears not to be restored at T c . We mention that the effective breaking 
of the axial U{\)a symmetry appears substantially reduced especially above T c , as inferred from 
the difference between the correlators in the pion and S channels [[16, 13 , [l4| p"5| . 



4 Again, this has been argued on the basis of a first-order e-expansion perturbative calculation 
in the framework of the corresponding effective three-dimensional quantum field theory. This has 
been recently confirmed by a six-loop computation within a three-dimensional perturbative scheme 

in. 
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bigous identification of the universality class of the finite-temperature transition in 
two-flavor QCD. Moreover, assuming that two-flavor QCD undergoes a continuous 
transition belonging to the three-dimensional 0(4) universality class, the 0(4) crit- 
ical equation of state gives the asymptotic relations, for T — > T c and in the limit of 
vanishing quark masses, among the reduced temperature t = T/T c — 1, the quark 
masses, and the expectation value of the bilinear Mj = (q^qRj). In particular, in 
the two-flavor case the relevant order parameter can be written as a four-component 
vector 0, whose expectation value M is related to the quark condensate (qq), and 
which is coupled to a uniform external field H related to the quark masses [|], ||, |3|, |] . 
The scaling properties can be written as 

M = B c H\H\ (1 - S)/s E(y), (1.2) 
y = {B/B c ) 1 ^t\H\- 1 '^\ 

where f3, 7, and 5 are the critical exponents, 5 B and B c are nonuniversal constants 
fixing the normalizations (they are the amplitudes of the magnetization at the coex- 
istence curve and along the critical isotherm, respectively), and E(y) is a universal 
function. Similar asymptotic relations can also be written for the (quark-mass) 
susceptibilities; they can be easily derived by taking appropriate derivatives of the 



relation (|1.2|) . Corrections to the scaling behavior (|1.2| ) are suppressed by powers 
of t and \H\. The leading ones are of order t A on the axis H = and of order 
|iJ| A /(/3+7) on the critical isotherm [in general corrections have the scaling behavior 



t A g(y)), where A « 0.58 |24|, |25] 



The critical equation of state of O(N) models has already been studied in e = 
4 — d expansion at two loops |27| and in 1/N expansion at order 1/N |[28|| . However, 
these results do not allow us to obtain a quantitatively accurate determination of the 
critical equation of state, essentially because the available expansions are too short. 
On the other hand, on the numerical side, rather accurate results have been obtained 
from Monte Carlo simulations of the 0(4) spin model ([Li] ) p9| , |30| , pT| . 

In this paper we compute the critical equation of state by using a different 
analytic method. The starting point is the small-field expansion of the effective 
potential (in statistical-mechanics language it corresponds to the small-magnetization 
expansion of the Helmholtz free energy), or equivalently of the equation of state, in 
the symmetric (high-temperature) phase (t > 0). The first few nontrivial terms of 
this expansion have already been determined by exploiting field-theoretical methods 



5 The critical exponents of the 0(4) universality class have been determined rather accurately by 
exploiting field-theoretical methods and lattice techniques. We mention the field-theory estimates of 
Ref. ||, v = 0.741(6) and r) = 0.035(5), and the Monte Carlo results of Ref. |§, v = 0.749(2) and 
7/ = 0.0365(10). A more complete list of references and results can be found in Ref. f2q| . The other 
exponents can be obtained by using scaling and hyperscaling relations: a = 2 — 7 = v(2 — 77), 
(3 = v{\ + rj)/2 and S = (5 — 77) /(l + 77). For example, using the most precise Monte Carlo results 
one obtains a = -0.247(6), 7 = 1.471(4), (3 = 0.388(1), and 5 = 4.789(6). 
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in the continuum <p theory [^, 33]. We use them to determine approximations of 
the equation of state that are valid in the whole critical regime. This requires an 
analytic continuation from the high-temperature phase (t > 0) to the the coexistence 



curve (t < 0), which can be achieved by using parametric representations [34, 35 



36], |37| implementing in a rather simple way the known analytic properties of the 
equation of state (Griffiths' analyticity) . We construct systematic approximation 
schemes based on polynomial parametric representations of the critical equation of 
state that have the correct analytic properties, take into account the Goldstone 
singularities at the coexistence curve, and match the known small-field behavior 
of the effective potential. Using the results for the critical equation of state, we 
then derive estimates of several universal amplitude ratios. As we shall see, we 
will provide rather accurate determinations of the scaling functions related to the 
magnetization, the longitudinal (quark-mass) susceptibility, and of several universal 
amplitude ratios, involving also quantitites defined at the so-called pseudo-critical 
line that corresponds to the maximum of the longitudinal susceptibility at fixed H. 

The paper is organized as follows. In Sec. |2| we discuss the general properties of 
the critical equation of state. In particular, we consider the small-field expansion of 
the effective potential and of the critical equation of state in the symmetric phase. 
In Sec. |] we describe our approximation schemes based on polynomial parametric 
representations, and determine the scaling functions that give the scaling behavior of 
the magnetization (quark condensate) and of the longitudinal susceptibility. Finally, 
in Sec. 01 we determine several universal amplitude ratios. 



2. The critical equation of state 
2.1 General properties 

—* — # 

The equation of state relates the magnetization M, the magnetic field H, and the 
reduced temperature t = (T — T c )/T c . In the neighborhood of the critical point t = 0, 
if = 0, it is usually written in the scaling form 

H = B~ 5 MM 5 ~ X f{x), x = t(M/B)-V p , (2.1) 

where M = \M\, f(x) is a universal scaling function normalized in such a way that 
/(— 1) = 0, /(0) = 1, and B c and B are the amplitudes of the magnetization on the 
critical isotherm and at the coexistence curve: 

M = B C H 1/S , t = 0, (2.2) 
M = B{-tY, t<0, H^0. (2.3) 

The scaling function f(x) and E(y), cf. Eq. (|L2|), are clearly related: 

E{y) = f{x)- l l\ y = xf(x)- 1 ^. (2.4) 
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The equation of state is analytic for \H\ > 0, and therefore f(x) is regular 
everywhere for x > — 1. In particular, f(x) has a regular expansion in powers of x 
around x = 



/(x) = l + ^/ n V, 

n=l 

and a large- X expansion of the form 

oo 

f(x)=X^fn^ 2nP - 



(2.5) 



(2.6) 



n=0 



At the coexistence curve, i.e. for x — * —1, the Goldstone singularities appear. 
General arguments predict that, at the coexistence curve and in three dimensions, 
transverse and longitudinal susceptibilities behave respectively as 



XT 



M 
Jf 1 



XL 



dM 
~dH 



[2.7) 



In particular, the singularity of \l for t < and H — > is governed by the zero- 



temperature infrared-stable Gaussian fixed point p8j, |38], |39fl , leading to the predic- 
tion 



/(x)^/ 2 coex (l+^ 



for 



x 



-1. 



:2.8i 



The nature of the corrections to the behavior (|2.8|) is less clear. It has been conjec- 
tured [flO, [|T], using essentially e-expansion arguments, that, for x — > — 1, i.e., 
near the coexistence curve, v = 1+x has a double expansion in powers of w = HM~ 5 
and w/ d ~ 2 )/ 2 . This would imply that in three dimensions /(x) could be expanded 
in integer powers of v at the coexistence curve. On the other hand, an explicit 



calculation J|2] to next-to-leading order in the 1/N expansion shows the presence 
of logarithms in the asymptotic expansion of f(x) for x — > —1, so that f(x) can- 
not be expanded in powers of v. However, these nonanalytic terms represent small 
corrections of order 0(t> 2 logt>) compared to the leading behavior ( |2.8j ). 

Finally, beside the scaling functions E(y) and f(x), we introduce a scaling func- 
tion associated with the longitudinal susceptibility, by writing 



XL 



where 



D(y) 



E(y) - ^E>(y) 



.1-1/6 



(2.9) 



(2.10) 



The function D(y) has a maximum at y — y r 
pseudocritical line t ma _ x (H), see also Sec. [|. 



(35f(x)-xf'(xY 
ax corresponding to the crossover or 
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2.2 The small-field expansion of the effective potential 

Let us consider the continuum yj 4 theory that describes the three-dimensional 0(4) 
universality class, i.e. 

S = J d z x\^d^{x)-d^{x) + K0{x)-0{x) + ^u[${x) ■ ${x)f - H ■ ${x)} , (2.11) 

where <p(x) is a four- component real field. Let us also consider a zero-momentum 
renormalization scheme [131 (see also Refs. |34], H[): 

r2»(p) = ^Z- 1 [m 2 +p 2 + 0(p 4 )] , (2.12) 

rSL(O) = Z ^ m \ (SabScd + SacSbd + S ad S bc ) , (2.13) 

where ri",... )0n are n-point one-particle irreducible correlation functions, m and g are 
the zero-momentum mass scale and quartic coupling, respectively. 

The effective potential (Helmholtz free energy) is related to the (Gibbs) free 
energy of the model. If M = (0) one defines 

JF(M) =M-H- -logZ(H), (2.14) 

where Z(H) is the partition function and the dependence on the temperature is 
always understood in the notation. The effective potential J-'(M) is the generator 
of the zero-momentum one-particle irreducible correlation functions. In the high- 
temperature phase it admits an expansion around M = 0: 

oo 

A jr = HM) _ m = £ -_„ 2; ,U^ (2.15) 

This expansion can be rewritten in terms of the expectation value of the renormalized 
field (p r {x) = Z^^ipfa), i.e. M r = (<p r ), 

AT = im 2 M 2 + m3 ~ j J^y92jM? j . (2.16) 

Note that g± is the renormalized coupling g appearing in Eq. ( |2.13| ); its critical 
limit is the fixed point of the renormalization-group equations, that is the zero of 
the Callan-Symanzik /3-function (3(g) = mdg/dm. The coefficients g^j approach 
universal constants (which we indicate with the same symbol) for m — > 0. By 
performing a further rescaling 

M 2 r =-z 2 (2.17) 
9a 

in Eq. (|2.16 ), the effective potential can be written as 

AF= — A(z), (2.18) 

94 
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where 



and 



4! 'j^m 

r 2j = 3 > 3- (2.20) 

The universal quantities #4 and the first few r 2 j have been estimated by using 
field-theoretical methods. Fixed-dimension perturbative calculations provide the 
estimates g± = 17.30(6) (from an analysis |M of the six-loop expansion of [3(g) 



), 7*6 = 1.81(3) (from an analysis p2]| of the corresponding four-loop series f33fl 



and r 8 = 0.456 (three-loop series |3~3fl), while e-expansion computations give |32] 
g4 = 17.5(3) (four-loop series), r 6 = 1.780(8) and r 8 = 0.2(4) (three-loop series). 
Other substantially less precise results and references can be found in Refs. f3"2" |. 



Since z oc t where t = T/T c — 1 oc r — r c , the equation of state can be 
written in the form 

dT(M) ^ M „ s 



with 



^)4- + f 3+ 5^ 2 "' (2 - 22) 

The two functions f(x) and F(z) are related: 

z~ s F(z) = F™f(x), z = zox'P, (2.23) 

where z = (R^) 1 / 2 is a universal constant, see Sec. fl[ Because of Griffiths' analytic- 
ity, J-(M) has also a regular expansion in powers of t for \M\ fixed. Therefore, F(z) 
has the large- z expansion 

F(z) = z 5 J2 F " z ~ k/ ^ (2-24) 

fc=0 

3. Parametric representations of the equation of state 
3.1 Polynomial approximation schemes 

In order to obtain a representation of the equation of state that is valid in the 
whole critical region, we need to extend analytically the expansion ( 2.19|) to the low- 



temperature region t < 0. For this purpose, we use parametric representations that 
implement the expected scaling and analytic properties. We set [KB 



M = m o R m(9), 
t = R(l-6 2 ), 
H = hoR ps h{B), (3.1) 
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-1 -0.5 0.5 1 

X 

Figure 1: The scaling function f{x) as obtained by the n=0, n=l (A), and n=l (B) 
approximations. For comparison, we also show the Monte Carlo (MC) result of Ref. . 
although it can be hardly distinguished from the n=l (A) curve. 



where h and mo are normalization constants. The variable R is nonnegative and 
measures the distance from the critical point in the (t, H) plane, while the variable 
9 parametrizes the displacement along the lines of constant R. The functions m(9) 
and h(9) are odd and regular at 9 = and at 9 = 1. The constants m and ho can 
be chosen so that m{9) = 9 + 0(9 3 ) and h{9) = 9 + 0(9 3 ). The smallest positive zero 
6*0 of h(9), which should satisfy 9 > 1, corresponds to the coexistence curve, i.e., 
to T < T c and H — > 0. The parametric representation satisfies the requirements of 
regularity of the equation of state. Singularities can only appear at the coexistence 
curve (due, for example, to the logarithms discussed in Ref. P%|), i.e., for 9 = 9 . 
The mapping (|3.1|) is not invertible when its Jacobian vanishes, which occurs when 

Y(9) = (1 - 9 2 )m\9) + 2(39m{9) = 0. (3.2) 



Thus, parametric representations based on the mapping ( |3.1| ) are acceptable only if 
9 < 9i, where 9\ is the smallest positive zero of the function Y{9). 

The functions m(9) and h(9) are related to the scaling function f(x) through 



1-9 



2 r 



x 



9 2 



m(flo) 
m(0) 



(3.3) 



/(*) 



m{9)y b h(9) 
m(l)J /i(l) 
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Figure 2: The scaling function F(z). 



The asymptotic behavior Q2.8|) is reproduced by requiring that 

n n\2 



h{0) 



for 9 -> 6 n . 



(3.4) 



Following Ref. P?| , we construct approximate polynomial parametric representations 
that have the expected singular behavior at the coexistence curve (Goldstone singu- 
larity) and match the known terms of the small-;? expansion of F(z), cf. Eq. ( |2.22| ). 
We consider two distinct approximation schemes. In the first one, which we de- 
note by (A), h{9) is a fifth-order polynomial with a double zero at 9q and m(9) is a 
polynomial of order (1 + 2n): 



scheme (A) : m(9) = 9 ( 1 + ^ gO 

i=\ 

h{e) = e(i-e*/el) 2 . 

In the second scheme, denoted by (B), we set 
scheme (B) : m(8) = 6, 

n 

h(e) = e(i-e 2 /e 2 o ) 2 (i + J2ci0 



(3.5) 



(3.6) 



i=i 



Here h(9) is a polynomial of order 5 + In with a double zero at 9q. For n=0 ap- 
proximations (A) and (B) coincide. It is worth noting that the n=0 approximation 
becomes exact in O(N) models for N oo [BR In both schemes the parameter 
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Figure 3: The scaling function E(y). 



8 and the n coefficients Cj are determined by matching the small- z expansion of F(z). 
The scaling function F(z) can be written in terms of the parametric representation 

as 

z = pm(6) (1 -9 2 )'\ (3.7) 
F(z(6))=p(l-6^ S h(6), 

where p is a normalization [[35], p6[| , which is fixed by matching the expansion ( |2.22|) 
to order z 3 . Explicitly, we have 

12 

p 2 = 6/5(5-1) - — -66ci, (3.8) 

where 6 = 1, — 1 respectively for scheme (A) and (B). In both schemes, in order to 
fix #o and the n coefficients q we need to know the values of n + 1 coefficients r2j, 
i.e., r 6 , ...r 6+2n - This method has been already applied to the determine the critical 
equation of state of the XY |J7|, [ITJ and Heisenberg [HS| universality classes. 



3.2 Results 

In order to implement the above-presented approximation schemes, we use the Monte 
Carlo estimates |25f v = 0.749(2) and t] = 0.0365(10) for the critical exponents, and 



the field-theoretical estimates r$ = 1.79(2) and r§ = 0.2(4) which take into account 
both fixed- dimension and e-expansion results, see Sec. |2.2| . This will provide three 
different approximations: n=0, n=l (A), and n=l (B). We find 

6l = 2.795(40) n = 0, (3.9) 
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Figure 4: The scaling function D(y). 



and 

6l = 2.949(150), ci = -0.0225(200), n = 1 (A), (3.10) 
9l = 2.373(200), ci = 0.0650(300), n = 1 (B), (3.11) 

where the number between parentheses indicates how much the coefficients vary when 
the input parameters change by one error bar. The relatively small value of c\ in 
the n=l approximations supports the effectiveness of the approximation schemes. In 
Figs. [I], ||, §, and H we show respectively the scaling functions f(x), F(z), E(y), D(y), 
as obtained from the n=0, n=l (A), and n=l (B) approximations for the central 
values of the input parameters. The differences among the three approximations 
give an indication of the size of the systematic error of the approximation schemes. 
The figures show that it is rather small, suggesting that the n — 0,1 representations 
already provide good approximations of the critical equation of state. 

In Fig. [I] we compare our results with the scaling function f(x) obtained in 
Ref. [[ItJ by interpolating Monte Carlo data for the 0(4) spin model ( |1 . 1| ) . The 



agreement is overall satisfactory. The Monte Carlo determination of f(x) turns out 
to be hardly distinguishable from the n=l (A) curve. A more direct comparison with 
the Monte Carlo data of Ref. is presented in Fig. |5|, where we plot Mj H 1 ^ versus 
(J-J c )/H l /P & using 5 = 4.789(6), (3 = 0.388(1) and the raw data reported there. To 
avoid finite-size and scaling corrections we only consider the results for the largest 
lattices that are closest to the critical point, in particular those satisfying H < 0.01 
and | J — J c \ < 0.05. The full line shown in Fig. ^] is the n=l (A) approximation, cf. 
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n = n = 1 (A) n = 1 (B) final estimates 



/o 


1.44(1) 


1.53(9) 


1.48(3) 


1.50(8) 


fa 


0.29(1) 


0.35(6) 


0.32(2) 


0.33(5) 


fa 


— 0.060(4) 


—0.08(1) 


a at/ i \ 

—0.07(1) 


—0.08(2) 


fo° 


0.82(1) 


0.92(10) 


0.87(3) 


0.89(9) 


fcocx 
J2 


4.8(3) 


2.2(1.3) 


3.3(6) 


2.8(1.4) 




-0.3(1) 


*0.2(4) 


*0.2(4) 




rio 


2.7(5) 


-4(4) 


-6(7) 


-5(6) 


rp co 
^0 


0.0236(3) 


0.0241(5) 


0.0240(4) 


0.0240(5) 


1/max 


1.46(1) 


1.35(10) 


1.40(4) 


1.38(9) 


D(y max ) 


0.3417(4) 


0.3429(11) 


0.3426(7) 


0.3427(10) 



Table 1: Results concerning the scaling functions obtained by using the n=0, n=l (A), and 
n=l (B) approximations. Numbers marked with an asterisk are inputs, not predictions. 



Eqs. ( |3~5| ) and fl3.10| ), while the dotted lines show the uncertainty due to the error on 



the corresponding input parameters. The necessary normalization conditions related 
to the amplitudes B and B c , cf. Eq. ( [L.2|) , have been determined by fitting the data. 



The agreement is good: most Monte Carlo data are within the uncertainty of the 
n=l (A) curve. 

In Table [l] we report results concerning the behavior of the scaling functions 
f(x) and F(z) for H = and on the critical isotherm, cf. Eqs. ( |2.5| ), ( |2.6|) , ( |2.8| ), 
( p,19| ), and ( |2.24| ). We also report estimates of y max and D(y max ), where y max is the 
value of y corresponding to the maximum of the scaling function D(y) defined in 
Eq. ( |2.9| ). The displayed errors refer only to the uncertainty of the input parameters 
and do not include the systematic error of the procedure, which may be determined 
by comparing the results of the various approximations. Our final estimates are 
reported in the Table [3]. They are obtained by taking the average of the results 
n—1 (A) and n—1 (B). The error is just indicative: it is the sum of the uncertainty 
induced by the input parameters (we take the average of the two errors) and of half 
of the difference between the two approximations. In most cases this procedure leads 
to an estimate that includes the n=0 result within its error. For comparison, we also 
report the estimates / 2 coex = 2.20(5), / °° = 0.888(7), F °° = 0.018(4), y max = 1.33(5), 
and -D(y max ) = 0.341(1), which can be obtained from the results of the fits of the 
Monte Carlo data reported in Refs. p0| , |3~l|j. 6 Finally, we mention that the result for 
rio is consistent with the estimate rio = 9(17) obtained from the analysis of its 0(e 3 ) 
series 



32 



6 The comparison of our results with the Monte Carlo estimates of Refs. |30| |T) should be done 
with caution. Monte Carlo results are subject to scaling corrections and finite-size effects. Moreover, 



Ref. (pfj used S — 4.86 and (3 — 0.38, while we use the more precise estimates 25 8 — 4.789(6) and 



P = 0.388(1). Our errors take into account the uncertainty on the critical exponents. 
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Figure 5: Monte Carlo data for M/H 1 / 5 versus (J — J c )/H l ^ & from Ref. |3(| compared 
with the n=l (A) approximation (full line). The dotted lines show the uncertainty on the 
curve due to the errors on the input parameters. The inset shows the differences between 
the Monte Carlo data and the n=l (A) curve: the errors are related to the uncertainty on 
the theoretical curve; the statistical errors of the Monte Carlo data are much smaller. 



Universal Amplitude Ratios 



U = A+/A- 
Rt = -CtB 2 /{C + f 
g, = -Q + /[(C + ) 2 (D 3 ] 
Rt = (Q + ) 1/3 

P m = T^B/B C 



Rt = aA+C+/B 2 
R x = C + B & - l /B 5 c 

Q c = B\rf/C + 
Rp = C + /C p 



Table 2: Definitions of several universal amplitude ratios. 



4. Universal amplitude ratios 



Universal amplitude ratios characterize the critical behavior of thermodynamic quan- 
tities that do not depend on the normalizations of the external (magnetic) field, of 
the order parameter (magnetization), and of the temperature. From the scaling func- 
tion f(x) one may derive many universal amplitude ratios involving zero-momentum 
quantities, such as the specific heat, the magnetic susceptibility, etc.... For example, 
the universal ratio of the specific-heat amplitudes in the two phases can be written 
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as 

^ = (41) 

where, in the case of the 0(4) universality class for which — 1 < a < 0, 

<p{x) = XNQ ~ 2 { /(0) + - \xr 2 f(x) + [ X dy \yr> [f(y) - f'(0) - yf"(0)} . 

a - 1 a J 

(4.2) 

Several universal amplitude ratios can be expressed in terms of the amplitudes de- 
rived from the singular behavior of the specific heat 

C H = A ± \t\- a + b, (4.3) 

where b is a nonuniversal constant, the magnetic susceptibility in the high-temperature 
phase 

x = NC + t~^ (4.4) 

with N = 4, the zero-momentum four-point connected correlation function in the 
high-temperature phase 



xt = m±Ac t t-^ (4.5) 

(again with N — 4), the second- moment correlation length in the high-temperature 
phase, corresponding to the inverse mass in the zero-momentum renormalization 
scheme, cf. Eqs. (|2.12|) and (|2.13|) , 

£ = f + t~\ (4.6) 



the spontaneous magnetization on the coexistence curve, cf. Eqs. ( |2.2| ) and ( |2.3| ). 
We also consider amplitudes along the crossover line t ma . x (H), that is defined as 
the reduced temperature where the longitudinal magnetic susceptibility XL{t,H) = 
dM/dH has a maximum at H fixed, i.e. 

t m UH) = T p H 1 ^\ (4.7) 
XL(t m ^H)=C p t m l x . (4.8) 

In Table ^| we report the definitions of several universal amplitude ratios that have 
been considered in the literature. Their expressions in terms of the functions m(8) 
and h{0) can be found in Ref. j48fl . Note the following relations 



fs° = r~ x \ /r = (4.9) 

X 

2/max P-rrl^ i D(jJ max ) R p P, m R^- (4.10) 



where f£° and /f° are related to the large-x behavior of f(x), cf. Eq. ( |2.6| ), and y n 
is the value of y where the scaling function D(y) takes its maximum. 
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n = n = 1 (A) n = 1 (B) final estimates 





c\ no / o\ 

2.03(3) 


1.89(11) 


1.93(6) 


1.91(10) 


R i 


0.250(5) 


0.274(25) 


0.263(9) 


0.27(2) 




8.04(7) 


7.5(5) 


7.7(2) 


7.6(4) 


R x 


1.22(2) 


1.09(12) 


1.15(4) 


1.12(11) 


p 


1.159(3) 


1.12(3) 


1.140(11) 


1.13(2) 


R p 


2.049(3) 


2.041(7) 


2.042(6) 


2.042(7) 




2.01(4) 


2.05(5) 


2.03(3) 


2.04(5) 



Table 3: Estimates of several universal amplitude ratios. 



In Table |3] we report the estimates of several universal amplitude ratios, as de- 
rived by using the approximations n=0, n—1 (A), and n=l (B). Again, the displayed 
errors refer only to the uncertainty of the input parameters. Our final estimates are 
determined as in the case of the quantities reported in Table [I]. They can be again 
compared with the Monte Carlo results reported in Refs. |30|, [JT|: R x = 1.126(9), 
Rf = 8.6(9) and P m = 1.11(2). There is good agreement, keeping again into account 
that Refs. |30|, |3lJ used different values of the critical exponents (see footnote ||). 



Moreover, using Eq. (|4. 1| ) and the approximate interpolation formula of Refs. [p0| , pl| 
we obtain Uq = 2.0(2), where the error is obtained by varying the curve parameters 
within their error. We also mention the estimates ||9|| Uq ~ 2.044 and R£ ~ 0.263, 



obtained by an analysis of three-loop series in the minimal-subtraction scheme with- 
out e expansion. They are in good agreement with those presented in Table ^. 

Universal amplitude ratios involving the correlation-length amplitude / + , such as 
R^ and Q c , can be obtained by using the estimate of g 4 . If we take = 17.30(6), 

we obtain 

(r>+ E>+ \ X / 3 
= 0.490(4), (4.11) 

Q c = — = 0.44(2). (4.12) 
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